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The genus Elaphe is represented by the following eleven 
taxa in Eurasia (Schulz 2013): E. anomala (Boulenger, 
1916), E. bimaculata Schmidt, 1925, E. carinata (Gunther, 
1864), E. climacophora (Born, 1826), E. davidi (Sauvage, 
1884), E. dione (Pallas, 1773), E. quadrivirgata (Boie, 
1826), E. quatuorlineata (Lacepede, 1789), E. sauromates 
(Pallas, 1811), E. schrenckii (Strauch, 1873), and E. zoige- 
ensis Huang, Ding, Burbrink, Yang, Huang, Ling, Chen 
& Zhang, 2012. 

While most Elaphe species in China can be easily distin¬ 
guished from each other, the Dione’s ratsnake, E. dione, can 
be confused with the twin-spotted ratsnake, E. bimaculata, 
because they closely resemble each other in their colora¬ 
tion and scalation. However, E. dione and E. bimaculata are 
parametric and treated as valid species due to differenc¬ 
es in colour pattern (Schulz 1996) and genetic distances 
based on 12s and coi gene sequences (Utiger et al. 2002). 

Dione’s ratsnake is the most widespread snake of the ge¬ 
nus Elaphe, occurring in Afghanistan, Armenia, Azerbai¬ 
jan, China, Georgia, Iran, Kazakhstan, the Kyrgyz Republic, 
Mongolia, North and South Korea, Russia, Tadzhikistan, 
Turkmenistan, Turkey, Ukraine, and Uzbekistan (Schulz 
2013). Even though E. dione has adapted to steppe habitats, 
it has a broad ecological range that permits it to inhabit a 
variety of habitats (Schulz 1996). In China, the species has 
been recorded from the following provinces and munici¬ 
palities: Anhui, Beijing, Gansu, Hebei, Heilongjiang, Hen¬ 
an, Hubei, Jiangsu, Jilin, Liaoning, Nei Mongol, Ningxia, 
Qinghai, Shaanxi, Shandong, Shanxi, Tianjin, and Xinjiang 
(Lu et al. 1999, Macey et al. 1988, Schulz 1996, Wallach 
et al. 2014). Reports from southern China, e.g., from Kun- 
shan and Shanghai, might be based on misidentifications 
of E. bimaculata (Schulz 1996). The species has been also 


reported from Sichuan, e.g., from the Tangjiahe National 
Nature Reserve (Chen et al. 1999), located at the north¬ 
eastern edge of Sichuan, from the Yalong and Dadu River 
(Schulz 1996), both flowing through the western part of 
Sichuan and from the Yele Nature Reserve at the south¬ 
western edge of the Sichuan Basin (Hu 2005) (Fig. 1). 

Here we report a new record of Dione’s ratsnake from a 
high-altitude area of Garze Tibetan Autonomous Prefec¬ 
ture, Sichuan, and provide genetic data on species identi¬ 
ty. We found two adult snakes and a complete exuvia at 
3,520 m a.s.l. (3i°4i’ N 99°44’ E, WGS84; Fig. 1A) in the 
north of the Chola Shan, near the Yalong River Valley on 
7 June 2014. The snakes were basking in the afternoon be¬ 
tween stones in a ruderal, vegetated area close to a small 
Tibetan village (Fig. 1B). No voucher specimens were col¬ 
lected. The dry exuvia was preserved in a standard 50 ml 
plastic centrifuge tube. 

Basic measurements of length and weight were record¬ 
ed in the field using a measuring tape and a spring balance 
to the nearest 1 cm and 1 g, respectively. Abbreviations are 
as follows: SVL (snout-vent length, measured from the tip 
of the snout to the anterior margin of the cloaca); TaL (Tail 
length, measured from the posterior tip of the anal plate 
to the posterior terminus of the tail button); and TL (total 
length, SVL+TaL). In addition, several scale counts were re¬ 
corded following Schulz (1996): VEN (ventrals); SC (sub- 
caudals); DSR (number of dorsal scale rows [1/2/3], counted 
at the level of one head length behind the head, at midbody 
[i.e., at the level of the ventral plate corresponding to half the 
total number of ventrals], and at one head length before the 
vent); IL (infralabials); SL (supralabials); P (postoculars). Bi¬ 
lateral scale counts were given as left/right. The two snakes 
were also buccal-swabbed using a PG-100 Performagene™ 
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kit in order to validate the species taxonomic identity by 
means of DNA-barcoding. DNA samples were extracted 
from the swabs with the PG-AC4 Performagene™ reagent 
package according to the manufacturer’s instructions. For 
DNA extraction from the skin sample, we used a peqGOLD 
MicroSpin Tissue DNA Kit (Peqlab Biotechnologie GmbH, 
Germany). DNA quantity and quality was assessed using the 
Epoch Microplate Spectrophotometer (BioTek, USA). DNA 
concentration in the buccal extracts was 23.93 and 38.28 ng/ 
pi, respectively, and 7.63 ng/pl in the skin extract. These val¬ 
ues are within the lower range of DNA concentrations usu¬ 
ally obtained from tissue samples, but similar or even higher 
than concentrations obtained using other non-invasive ap¬ 
proaches in vertebrates (e.g., faeces or hairs; Ball et al. 2007, 
Morin et al. 2001). The three DNA samples have been de¬ 
posited in the Natural History Museum Erfurt (SH1406T0, 
SH1406T1, SH1406T2). 

We obtained partial sequences of the mitochondrial 
12s rRNA subunit, nadh dehydrogenase subunit 4 (nd4), 
cytochrome b (cyt-b), and cytochrome oxidase subunit 
1 (coi) genes and the nuclear oocyte maturation factor 
Mos c-mos from the three samples. We used the primers 


L1091 and H1478 (Kocher et al. 1989) to amplify the 12s 
region, LEU and nd4 (Arevalo et al. 1994, Forstner et 
al. 1995) for the nd4 fragment, L14724 (Irwin et al. 1991) 
and Hi5547b (Kelly et al. 2003) for the cyt-b fragment, 
coidegi and coibdeg (Utiger et al. 2002) for coi, and S77 
and S78 (Lawson et al. 2005) for the c-mos region, respec¬ 
tively. Amplicons were purified using a QIAquick PCR Pu¬ 
rification Kit (Qiagen, Germany) and sequenced in both 
directions with the same amplification primers by Euro- 
fins-MWG-Operon (Ebersberg, Germany). Sequences of 
eight (nd4, cyt-b, c-mos) or eleven (12s, coi) Elaphe spe¬ 
cies were derived from GenBank (Supplementary table 1) 
and aligned with sequences from our samples using Bio- 
Edit 7.09.0 (Hall 1999). Protein-coding gene fragments 
(nd4, cyt-b, coi, c-mos) were translated into amino acids 
and no stop codons were observed. No heterozygous posi¬ 
tion was present in the c-mos sequences. To assess the evo¬ 
lutionary distance between these species and our samples 
we calculated uncorrected p-distances among haplotypes 
of the four partial genes in Mega 6.06 (Tamura et al. 2013), 
using the pairwise deletion option and considering both 
transitions and transversions. New sequences generated 
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Figure 1. Sampling site and habitat of Elaphe dione in Garze Prefecture, Sichuan. A) The record site of E. dione north of the Chola 
Shan is indicated by a star (map after Miehe et al. 2001). Green dots represent reports of E. dione according to Schulz (1996) and 
Hu (2005). The type locality and another record of E. zoigeensis have been also included (Huang et al. 2012, Schulz 2013). B) The 
white arrow points to the location where we found specimens of E. dione. Photo by S. Hofmann. 
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Table 1. Measurements and scale counts of E. dione from Sichuan. 



Adult individual 
(sample 
SH1406T1) 

Adult individual 
(sample 
SH1406T2) 

Exuvia 

(sample 

SH1406T0) 

Weight (g) 

184.0 

159.0 

- 

SVL (mm) 

740.0 

680.0 

850.0 

TaL (mm) 

140.0 

140.0 

200.0 

TL (mm) 

880.0 

820.0 

1050.0 

VEN 

not counted 

not counted 

188 

SC 

58 

58 

62 

DSR 

not counted 

not counted 

25-22-19 

IL 

10/10 

11/10 

not counted 

SL (eye contact) 

8/8 (2/2) 

9/8 (2/2) 

8/8 (2/2) 

P 

3/3 

2/2 

2/3 


in this study were deposited in GenBank (Supplementary 
table 1, accession codes KP091848, KP091851, KP091850, 
KP091853, KP115289, KP115290). 

Morphological characteristics and measurements are 
summarized in Table 1. These data, accompanied with fea¬ 
tures of coloration and pattern (Figs. 2A-E), fit in with 
those reported for E. dione (Schulz 1996). Since detailed 
descriptions of E. dione have been published frequently in 
the past, we omit the descriptions of the colouration and 
body/headshapes of the two individuals. 

Uncorrected p-distances ranged from -4.1-11.5% for 12s, 
10.9-15.5% for nd4, 9.8-13.7% for cyt-b, 8.2-15.0% for coi 
(Supplementary tables 2 and 3), and between 0.2% and 1.1% 
for c-mos (data not shown) between species. Intraspecif¬ 
ic divergence ranged from o.0-2.1% for 12s, 0.0-23% for 
coi, or was lower than 4.0% (nd4) and 3.0% (cyt-b), re¬ 
spectively (Supplementary tables 2 and 3). The close rela¬ 
tionship between E. dione and E. bimaculata is highlight¬ 
ed by the generally low levels of genetic distance between 
them compared to the p-distances between other Elaphe 
species. However, the extremely low divergences between 
sequences of E. bimaculata that were derived from the 
complete mitochondrial genome data (GenBank accession 
codes NC_024743; Yan et al. 2014) and all other sequences 
of E. dione used in our study suggest that the former are 
based on sample misidentification and were probably de¬ 
rived from E. dione specimens too. There is also striking 
molecular evidence that E. anomala and E. schrenckii, for¬ 
merly classified as subspecies E. schrencki anomala (Bou- 
lenger, 1916), might be conspecific (Supplementary ta¬ 
ble 2). Moreover, the sequence divergences between the 
E. sauromates specimens are of similar magnitude com¬ 
pared to the interspecific distances between E. sauromates 
and E. quatuorlineata, emphasizing the need for further 
studies of the taxonomic status of this species complex 
(Supplementary table 2). Our results also support that the 
recently described E. zoigeensis represents a (phylogeneti- 
cally) distinct species. Although its morphological charac¬ 
ters seem very similar to those reported for E. dione, the 
genetic distances between E. zoigeensis and all other Elaphe 


species warrant its species status under an evolutionary 
species concept (de Queiroz 2007). The low levels of 12s 
rRNA gene sequence divergence between our sample and 
the E. bimaculata and even the E. carinata specimens (Sup¬ 
plementary table 2: no. [2] to [4] and [7], respectively) are 
difficult to interpret, in particular because these specimens 
were purchased through the pet trade without locality 
data. In contrast, with respect to the coi data, the sequence 
distances between these specimens and our samples were 
within the range of interspecific divergence. Interestingly, 
we observed higher p-distance values between sequences 
of E. dione and our samples for all mt-markers and par¬ 
ticularly for coi (5.1-5.8%; Supplementary table 2) and nd4 
(3.7%; Supplementary table 3). Remarkably, the coi-differ- 
ences between E. dione samples from locations as far apart 
as Korea and Ukraine (Supplementary table 2 no. [12] and 
[14], respectively) are smaller compared to those observed 
between our samples and the Ukrainian sample. Overall, 
however, at least for coi, the p-distance values between 
E. dione and our samples are substantially lower compared 


A 



Figure 2. Head of Elaphe dione. A) Dorsal; B) frontal; C) right 
lateral; and D) left lateral views; E) back of the body. 
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to the mean distances between groups of species (8.2% 
to 13.8%; excluding the misidentified E. bimaculata and 
grouping E. anomala and E. schrenckii together; data not 
shown). Furthermore, we found no evidence of nDNA di¬ 
vergence between our sample and E. dione (c-mos, p-dis- 
tance 0.0%, data not shown). 

The area from which our samples originated belongs 
to the Hengduan Shan region that consists of a series of 
north-south aligned, highly rugged mountain chains al¬ 
ternating with deep gorges and river valleys, resulting in a 
complex topography, heterogeneous climate, and ecologi¬ 
cal niches. During the Quaternary glaciations, the Heng¬ 
duan Shan formed a mosaic of temporarily isolated ref¬ 
uges, contributing to the currently observed exceptional 
intra- and interspecific diversity and a high level of ende¬ 
mism across western Sichuan and northeastern Yunnan 
(Gao et al. 2007, Lei et al. 2007, Lopez-Pujol et al. 2011, 
Huang et al. 2013, Wen et al. 2014, Zhang et al. 2009). 
Today, that mountain area is recognized as a unique bio¬ 
diversity hotspot, emphasising its biogeographical signifi¬ 
cance (Boufford & van Dyck 1999, Myers et al. 2000, 
Boufford et al. 2004). It is reasonable to assume that the 
population from which we drew our samples is somewhat 
geographically isolated from other populations, which may 
result in regional molecular differences due to restricted 
gene flow. Moreover, since mtDNA is haploid, maternal¬ 
ly transmitted and, thus, has a smaller effective population 
size than does nuclear DNA, it is susceptible to stronger 
genetic drift, resulting in increased genetic divergence. The 
higher p-distances in mtDNA and the lack of differentia¬ 
tion in nDNA between our samples and E. dione are con¬ 
sistent with substantial intraspecific mtDNA sequence di¬ 
vergence due to genetic drift triggered by isolation and a 
potentially small population size, rather than (cryptic) spe¬ 
cies- or subspecies-level differentiation. 

Our sample site of E. dione in Garze confirms the pres¬ 
ence of the species in Sichuan, in the area north of the Cho- 
la Shan. There are a few previous reports of the species from 
areas around that site (reviewed in Schulz 1996, Hu 2005), 
indicating that the species is more common in Western Si¬ 
chuan than would be expected. Although E. dione has been 
reported to occur at altitudes of up to 3,500 m (Bannikov 
et al. 1977), our observation of the species at 3,520 m might 
be the first validated record for its occurrence at such high 
altitudes. The only known species from similar altitudes in 
Sichuan is Thermophis zhaoermi, which might be sympat- 
ric with E. dione at some locations. Given the ecological 
plasticity of the species, the presence of E. dione in areas 
farther south (e.g., in mountainous regions in northwest¬ 
ern Yunnan) is conceivable and should be investigated in 
the future. 
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Supplementary table 1 . Sequences for the genus Elaphe retrieved from GenBank and used in the present study. * Tissue collection of 
N. Helfenberger (Utiger et al. 2002); ** Yan et al. (2014); ‘The three samples yielded identical sequences for the respective partial 
genes except for col (variation for a single base pair at position 317 that could be A or G). 


Taxon 

Voucher/ isolate No. 

Origin 

12S 

GenBank accession codes 
nd4 cyt-b col 

c-mos 

Elaphe anomala 

(SH577*) 

China 

AY122803 



AY122719 


E. bimaculata 

(SHI 164*) 

(pet trade) 

AY122767 



AY122683 


E. bimaculata 

(SH1420*) 

(pet trade) 

AY122768 



AY122684 


E. bimaculata 

(SH1438*) 

(pet trade) 

AY122837 



AY122753 


E. bimaculata 

No voucher 

unknown 


DQ902283 

DQ902104 


DQ902062 

E. bimaculata 

(complete genome)** 

unknown 

NC_024743 

NC_024743 

NC_024743 

NC_024743 


E. carinata 

(SH1177*) 

(pet trade) 

AY122838 



AY122754 


E. carinata 

(SH575*) 

(pet trade) 

AY122839 



AY122755 


E. carinata 

LSUMZ 37012 

China 


DQ902284 

DQ902133 


DQ902063 

E. climacophora 

(SH528*) 

Japan 

AY122770 



AY122686 


E. climacophora 

CAS 163993 

Japan 


DQ902285 

DQ902105 


DQ902064 

E. davidi 

(SH576*) 

China 

AY122775 



AY122691 


E. dione 

LSUMZ 45799 

Russia 


DQ902287 

DQ902107 


DQ902066 

E. dione 

(SH2166*) 

Kazakhstan 

AY122832 



AY122748 


E. dione 

(SH1389*) 

Kazakhstan 

AY122834 



AY122750 


E. dione 

(SH1152*) 

Korea 

AY122829 



AY122745 


E. dione 

(SH1155*) 

Russia 

AY122830 



AY122746 


E. dione 

(SH545*) 

Ukraine 

AY122831 



AY122747 


E. dione 

(SH952*) 

China 

AY122833 



AY122749 


E. dione (3 samples)* 

No voucher, SH1406T0 

this study 

KP091848 

KP091851 

KP091850 

KP115289 

KP091853 

E. dione 

No voucher, SH1406T2 

this study 




KP115290 


E. quadrivirgata 

No voucher 

Japan 


DQ902300 

DQ902120 


DQ902078 

E. quadrivirgata 

(SH551*) 

Honshu, Japan 

AY122793 



AY122709 


E. quadrivirgata 

(SH470*) 

Honshu, Japan 

AY122794 



AY122710 


E. quatuorlineata 

LSUMZ 40626 

Turkey 


AY487067 

AY486931 


AY486955 

E. quatuorlineata 

(SH1112*) 

Italy 

AY122796 



AY122712 


E. quatuorlineata 

(SH975*) 

unknown 

AY122798 



AY122714 


E. sauromates 

(SH972*) 

Turkey 

AY122795 



AY122711 


E. sauromates 

(SH556*) 

Turkey 

AY122797 



AY122713 


E. schrenkii 

No voucher 

unknown 


DQ902302 

DQ902124 


DQ902082 

E. schrenkii 

(SH857*) 

(pet trade) 

AY122804 



AY122720 


E. zoigeensis 

isolate ZJ200801-02 

China 

HQ330522 

HQ330526 

HQ330524 



E. zoigeensis 

isolate ZJ200803 

China 

HQ330523 

HQ330527 

HQ330525 



E. zoigeensis 

HS2010014 

China 




JF510474 


E. zoigeensis 

HS2010016 

China 




JF510475 



I 
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[23] E. schrenckii AY122804/AY122720 0.3* 5.3 5.6 5.6 5.0 4.4 4.4 5.6 5.6 5.0 5.0 6.5 5.0 4.7 5.6 5.6 4.7 4.7 5.6 5.3 5.3 6.5 12.7 12.7 

[24] E. zoigeensis HQ330522/JF510474 7.7 8.0 8.0 7.7 9.1 7.9 8.2 9.1 10.9 8.5 8.5 9.4 9.1 9.4 9.1 7.6 9.1 9.1 9.7 9.4 10.0 11.5 8.0 0.4 

[25] E. zoigeensis HQ330523/JF510475 6.8 7.1 7.1 6.8 8.2 7.1 7.4 8.2 10.0 7.6 7.6 8.5 8.2 8.5 8.2 6.8 8.2 8.2 8.8 8.5 9.1 10.6 7.1 0.9 





































































Online Supplementary data - A new record of Elaphe dione from Western Sichuan 


Supplementary table 3. Uncorrected p-distances of the partial cyt-b gene (779 bp; lower left) and partial nd4 gene (682 bp; upper 
right) for 11 sequences from 8 Elaphe species; species name and GenBank accession code. The p-distance between our samples and 
E. dione is highlighted in bold. Conspecific distances are shadowed in grey. *[2] might be conspecific with E. dione [5] and our samples 
[6] or may indicate subspecies level; see also the annotations to Supplementary table 2. ** Higher distance values between E. dione [5] 
and our samples [6] indicate, again, differentiation due to genetic drift within the species’ range. 



[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 

[8] 

[9] 

[10] 

[11] 

[1] E. bimaculata DQ902104/DQ902283 


10.9 

12.5 

12.9 

11.1 

11.4 

13.6 

12.2 

14.1 

13.5 

13.6 

[2] E. bimaculata NC_024743 

10.0 


14.5 

13.3 

1.2* 

3.4* 

13.2 

13.7 

14.8 

14.2 

14.4 

[3] E. carinata DQ902133/DQ902284 

12.0 

11.3 


11.9 

14.2 

13.9 

11.0 

12.5 

12.3 

12.6 

12.8 

[4] E. climacophora DQ902105/DQ902285 

11.3 

11.4 

10.7 


13.0 

14.1 

12.2 

12.5 

12.9 

13.6 

13.8 

[5] E. dione DQ902107/DQ902287 

10.2 

0.4* 

11.4 

11.6 


3.7** 

12.9 

13.4 

14.8 

13.6 

13.8 

[6] E. dione (this study) KP091850/KP091851 

9.8 

2.6* 

12.1 

11.9 

2.4 


13.2 

13.1 

15.5 

15.0 

15.1 

[7] E. quadrivirgata DQ902120/DQ902300 

12.7 

11.3 

11.2 

10.8 

11.4 

11.8 


13.5 

12.8 

12.8 

12.9 

[8] E. quatuorlineata AY486931/AY487067 

11.8 

10.7 

12.5 

10.8 

10.8 

11.7 

12.3 


14.0 

12.8 

13.0 

[9] E. schrenckii DQ902124/DQ902302 

11.8 

11.7 

13.0 

11.3 

11.8 

11.4 

11.3 

11.7 


13.6 

13.8 

[10] E. zoigeensis HQ330524/HQ330526 

12.5 

10.9 

13.7 

11.6 

11.0 

11.4 

12.8 

11.3 

12.3 


0.11 

[11] E. zoigeensis HQ330525/HQ330527 

12.5 

10.9 

13.7 

11.6 

11.0 

11.4 

12.8 

11.3 

12.3 

0.0 
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